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Abstract 

Exploiting the "natural" frame of space curves, we formulate an intrinsic 
dynamics of twisted elastic filaments in viscous fluids. A pair of coupled 
nonlinear equations describing the temporal evolution of the filament's com- 
plex curvature and twist density embodies the dynamic interplay of twist and 
writhe. These are used to illustrate a novel nonlinear phenomenon: geometric 
untwisting of open filaments, whereby twisting strains relax through a tran- 
sient writhing instability without performing axial rotation. This may explain 
certain experimentally observed motions of fibers of the bacterium B. suhtilis 
[N.H. Mendelson, et al, J. Bacteriol. 177, 7060 (1995)]. 
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As D'arcy Thompson emphasized nature abounds with examples of the balance 
between growth and form. Thompson was especially interested in 'dynamic morphology,' 
the realization of this harmony in organisms. The balance is particularly manifest in the 
low Reynolds number regime of slow motions without inertia, where growth is generally 
dominated by dissipation and form is constrained by elasticity. The coupling between these 
two has been termed "elastohydrodynamics" . 

Possibly the most common 'form' of subcellular dynamics is the elastic filament: rod- 
like macromolecules perform a multitude of structural and biological roles. To gain some 
intuition for dynamics at this scale, the case of the elastic filament is thus a fruitful starting 
point 0. As a motivating example, we recall Mendelson's discovery of shape instabilities 
exhibited by a mutant strain of the bacterium B. subtilis Filaments of cells form as 
the rod-shaped bacteria fail to separate upon dividing; beyond a critical length, a filament 
buckles, twisting around itself to form a plectoneme (See Fig. |l|). It is thought that this 
instability arises from nonequilibrium twisting stresses generated in the cell walls through 
the process of growth 0,0. As the bacteria continue to grow, this process is iterated, leading 
to a hierarchy of supercoils. The handedness of the coils may be tuned with temperature or 
ionic concentrations. Since the filament ends are generally free and thermal fluctuations are 
negligible, this supercoiling is fundamentally different from that of overtwisted DNA that is 
closed or has pinned ends 0. As a first step toward understanding these highly nonlinear 
phenomena, we develop here an intrinsic deterministic formulation of twisted elastic filament 
motion, and apply it supercoiling of B. subtilis. Various elements of this formulation have 
appeared in previous work on closed elastic rods without dissipation 0], where a dynamic 
(rather than energetic [0) coupling between twist and bend degrees of freedom has been 
recognized 0. In the inertialess geometric argument for this coupling has been 

recently proposed and its qualitative consequences have been explored [jTU[. In presenting 



a general formulation of filament dynamics, we explicitly derive this twist-bend coupling and 
illustrate the qualitatively new physics of open filaments in a viscous medium. 

An important feature of our treatment is the parameterization of the filament. Although 
the "Euclidean" representation based on the position vector r(s) is convenient for linear 
stability analysis about simple shapes, it makes the full nonlinear problem unnecessarily 
complicated. As exploited by Kirchhoff the Euler angles of rigid body motion are 
better suited for elastic equilibria. This parameterization singles out a particular space- 
fixed reference frame and thus is partly extrinsic. An intrinsic representation expresses the 
rate of change of a frame {e.g. the material or Frenet-Serret frames) at each point of the rod 
in terms of the frame itself |T^, and is simplified by using the "natural frame" associated 



to the centerline. First used by Darboux |TB| and later by Hasimoto in his study of vortex 



filament motion 14 , this frame has been applied to the dynamics of space curves with 



bending but no twisting elasticity 



The configuration of an elastic rod may be described by a material orthonormal frame 



{ei, 62, 63} |jT2[: we choose 63 to point along the tangent, ei to point from the centerline to 
a painted stripe on the rod surface that connects one end to the other, and 62 = 63 x ei. 
The strains are given by the rate of rotation of this frame with arclength s, ds^i = flxei. 
To quadratic order in Q, the elastic energy of a homogeneous isotropic rod is 

£ = ^J ds{nl + nl) + ^J dsnj - J dsA, (i) 
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where Qi = Q ■ Bi, A and C are the bend and twist elastic constants respectively, and A(s) 
is a Lagrange multiplier function enforcing local inextensibility. 



A "natural" orthonormal frame |]T3| of a given curve consists of the unit tangent and 
two other vectors constrained to have zero tangential rotation rate. Imagine a weather vane 
sliding along the curve and always lying in the normal plane . If the weather vane spins 
around 63 with vanishing initial tangential angular momentum then it traces out one member 
of a family of natural frames related by s-independent rotations. Such a frame is constructed 
simply by untwisting the material frame, i.e., rotating ei and e2 around 63 by minus the 
accumulated twist angle -^{s) = f^ds'Q^. Since e* • c^^e = 0, the vector e = (ei + ^62) expi{} 
does not twist around 63. 

The rotation rate of this frame is given by the complex curvature ip: dsS = —ip^a 
dge^ = Iie{ipe*), where 

^p = {—Oil + ^22) expi J ds'^l3. (2) 

Some elementary shapes have simple ■^/'-representations: the straight line has = 0, the 
circle ip = a, and the helix ip = ae*^'^ (a, p are real constants). Likewise, certain geometric 
properties of curves are easily stated in terms of ip: a curve lies in a plane if and only if 
{Iieip{s),lmip{s)) lies on a line through the origin in the complex ip-plane, and a curve lies 
on a sphere if and only if {Iieip{s), lmip{s)) lies on a line not through the origin [|T7| . 

We compute from the elastic forces and moments per unit length using the principle 
of virtual work. If F is the force acting on a cross-section, the force per unit length Fs 
is found by varying the position of the centerline without rotating any element of the rod 
about the tangent direction: 6£ = — J dsFg ■ Sr{s), with the rotation angle 5% = 62 ■ 6ei = 0. 
Likewise, if M is the moment acting on a cross section, the tangential component of the 
moment per unit length Mg is found by varying the orientations of the elements of the rods 
without moving the centerline: 6S = — J dsM^ ■ e^Sx, with Sr = 0. Using ^r^ = {6r)s, 
Sea = {5x)<^abh - (ea ■ 5rs)e3 for a, 6 = 1, 2 and de^ = 6rs - (63 ■ 6rs)e3, we find 

dip = —2'ipe3 ■ 6rs + e ■ 5rss — iiplin. J ds'ipe ■ ^r^, (3) 
= 6xs - ^363 ■ + i^iBi + ^262) ■ Sts. (4) 

Note that keeping 6x = during the variation does not imply that the twist rate f2 = (^3 
stays fixed! (For brevity we drop the subscript on ^3 below.) 

Thus the force per unit length is —6S/6t = —e-^Ag + Re (ejFj_), where 

^± = -A{iPss + ^l^r^) + ^C{Q1P)^ - AiP, (5) 
and a term has been absorbed into A. Separating into its real and imaginary 



parts yields the classic results of Love [T^. Note that the force density is invariant under 
ip I— s- e^'^ip, e I— > e^'^e for constant (p. The tangential moment per unit length is —6S/6x{s) = 
Cfls- The conditions of zero force and moment at free ends imply the boundary conditions 

iP = ip^ = A = n = o. 

The dynamics of ip and fl follow immediately from (^) by choosing 5r = rtSt and Sx = 
Xt^t. Their form is completely dictated by geometry. (That for ipt is equivalent to earlier 
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results the undamped analogue of the Qt equation appears in The right-hand 

side of (§) expresses a fact easily demonstrated with the nearest rubber tube; there are two 
distinct ways of changing the twist density — (i) nonuniform axial rotations of the rod, and 
(ii) nonplanar bending motions. 

To complete the equations of motion, we use the simplest Rouse dynamics in which the 
only effect of the surrounding fluid is to provide a local isotropic drag ||19[ proportional to 
the local velocity: fdrag = C^"*- Likewise, we take the tangential moment per unit length to 
be mdrag = CrXt- Since at zero Reynolds number the elastic forces and moments must exactly 
balance the drag forces and moments, the complex curvature and twist obey 

CiJt = {d^s + + ^^Im fds'r^*^s - As^s (6) 

CQt = CQss + (Cr/C)Im^±.> • (7) 

Local inextensibility, dtds = 0, is assured by demanding A^s = Re(V^*jFj_). Eq. (^) is 
covariant and (|^) invariant under s-independent rotations of the natural basis about the 
tangent: ip i— >■ e^'^'ijj. Since the elastic energy (P is achiral, the dynamics are covariant under 
mirror reflections: ip ^ ip* , ^ —Vt. 

The derivation of the above dynamics relies solely on geometry and mechanics. However, 
for closed curves, we can reconsider (0) in light of the celebrated relation Ck = Tw + Wr 
||2CI|| . Here Ck is the linking number of the two curves traced out by the heads of the vectors 
ei and 62 (assuming a continuous material frame), Tw = J dsf2/(27r), and the writhe Wr 
is an integer plus A/{27i), where A is the signed area enclosed by the curve swept out on 
the unit sphere by 63 ||2T]|. For a locally inextensible filament, the twist dynamics (|^ can 
be interpreted as a local link conservation law under the assumption that smooth filament 
motions will not change the integer contribution to the writhe 0. The time derivative of 
the density wr of the noninteger contribution to the writhe is dtwr = ■ dte^ x dge^/ (27r) = 
— Im(jF^ *?/')/ (27r(^). Therefore, the twist dynamics has the form pt + js = 0, with the "link 
density" p = Q/{27t) + wr and link current j = — (C/Cr)^^- Of course, since the link is not 
a single integral of a local density there is no true link density; however, changes in link 
can be written as a single integral over a density [^]. While Ck and Wr are conventionally 
only defined for closed filaments, the equations of motion (|) and (0) show that the tradeoff 
between twist and writhe still has meaning for open filaments. 

With the coupled evolution equations for ijj and Q we turn to linear and nonlinear aspects 
of supercoiling instabilities. The simplicity of the ^/'-formulation for open filaments is seen 
first in the stability of a straight rod with a uniform twist density (temporarily neglecting 
the boundary conditions on Q). The linearization of (^) is 

(ipt = -Aipssss + iCQtpsss ■ (8) 

The eigenfunctions Wn of —Ad^+iCfldg have eigenvalues ct„ = —Ak^+CQk^ (and associated 
growth rates Un/O which are real since the operator is Hermitian. If = 0, the modes are 
planar shapes with Wn = c^^^ sin(fc„s) -|- c^^^ cos(/cns) -f- c^^^ sinh(/c„s) -|- c'^^ cosh(fc„s), with 
cos(fc„L) cosh(A;„L) = 1 Since the Wn{s) satisfy the boundary conditions and form 

a complete set, the shape of any free elastica in three dimensions, whether tightly wound 
helices or plectonemes, is expressible as 
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^(s,t) = £c„(t)iy„(s) . (9) 



n=l 



The time evolution is a nonlinear dynamical system in C„(t) = (Wnlipit)) , where the Wn 
are normalized. 

The linear stability analysis of a free rod with uniform twist is identical to that of 
a clamped rod with uniform twist. Using the complex Monge representation ^(z) = 
x{z) + iy{z), where r{z) = {x{z),y{z), z), yields C6 = -Mzzzz + iC^l^zzz with = = 
at the clamped ends. Note however that the Monge representation is incapable of describing 
plectonemes. The critical condition for the first centerline instability is Qc ~ 8.98A/CL [jl2|. 
When the twist is non-zero, each eigenfunction is a superposition of four helices with appro- 
priate spatial damping: Wn{s) = c^^-* exp{ik^\s), where for each n the four k!^^ are in 
general complex. Modes with negative eigenvalue ct„ have three helices of one handedness 
and one of the other; otherwise, all have the same handedness. Shapes calculated from the 
first two modes are shown in Fig. ^. 

The writhe-tracking term in (|^ is explicitly nonlinear; when the diffusion constant 
C/Cr = 0, twist only changes when the filament writhes. In assessing the relative im- 
portance of twist diffusion and writhe-tracking we note that while (r in (0) was interpreted 
as a rotational drag coefficient, in the context of bacterial filaments it is not clear that this 
dissipative mechanism dominates others, e.g., rearrangement of the polymer network inside 
the cell walls. If the shape instabilities are generated by nonequilibrium twist, then twist 
must not relax on time scales shorter than the observed buckling time (tens of minutes 
2^). Using the purely hydrodynamic twist diffusion constant C/Cr ||T0| , |23[| (with C = 10~^° 
erg-cm and Cr = 10^^° erg-sec/cm for a one micron diameter slender rod in water). 



and the typical filament length L ~ 100 fim for buckling, the time scale td = L'^(/C for 
twist to diffuse out the ends of the filament is only 10~^ sec, eliminating hydrodynamic 
drag as a candidate mechanism. This and other observations {e.g. of slow stress relaxation 
^) strongly suggest a dominant dissipative mechanism within the filament, characterized 



by a very much larger (r- Therefore, we consider the limit of negligible diffusion constant 
C/(r, where twist changes only by bending rather than by axial rotation of the filament, a 
relaxation mechanism we term "geometric untwisting." 

Geometric untwisting is illustrated in Fig. ^, obtained by numerical integration of 
(^) and (|^. An additional nonlocal contact force in !F± prevents self-crossing. The initial 
condition is a nearly straight filament with twist density of the form shown in Fig. ^a. 
This mesa-like profile reflects uniform local twist production during growth and diffusive 



boundary layers at the ends Motion begins as chiral buckling; as the curvature increases, 
twist is converted to writhe, and Q decreases in the interior (Fig. |^a). When twist has 
sufficiently relaxed and writhe is of order unity, the filament relaxes back to the straight 
ground state as link flows out the ends. In this overdamped dynamics the total elastic 
energy is monotonically decreasing, despite the transient increase in the bending energy 
during loop formation (Fig. ^d). 

The sequence in Fig. ^ is quite similar to that exhibited by short filaments of B. subtilis 
||2^ : after beginning the transition to a plectoneme via a looping instability, the ends touch 
only transiently, the loop snaps open, and the filament straightens. Plectoneme formation 
occurs after several aborted looping events, presumably as the twist stresses continually 
increase with filament extension. 
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The results presented here proceed from the robust to the specific: from kinematics 
dictated solely by geometry, through viscous dynamics governed by elasticity, and united in 
a description of the onset of bacterial supercoiling. This particular union will help unravel 
the puzzles surrounding growth-induced iterated writhing instabilities of bacterial fibers 
p6| and elucidate possible twist-induced instabilities in related systems, such as thermally 
fluctuating nucleic acids |27[| . 

We are indebted to Neil Mendelson and Charles Wolgemuth for ongoing collaborations, 
and thank A. Goriely, G. Huber, R. Kamien, J. Kessler, and M. Tabor for important dis- 
cussions. This work was supported by NSF Presidential Faculty Fellowship DMR96-96257 
(REG). 
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FIGURES 




FIG. 1. Supercoiled filament of B. subtilis [3]. Scale bar is 10 /im. Image courtesy of N.H. 
Mendelson. 




FIG. 2. Filament shapes at various twist densities. (a)-(c): ip = lOWi at 
= 0.22,1.0,2.22; (d)-(f): V = 8W2 at O/Q^^^ = 0.13,0.58,1.29. 
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FIG. 3. Geometric untwisting. Filament shapes are calculated from numerical integration of 
@ and (0) with C j A = 1. Time proceeds from left to right. 
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FIG. 4. Details of geometric untwisting, (a) Twist density evolution, with dotted line showing 
initial profile. Asymmetry in arises from that of the initial shape, (b) Time evolution of twist, 
bend, and total elastic energy. 
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